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SUMMARY 


Aerodynamic parameters of the High-Angle-of-Attack Research Aircraft 
(HARV) were estimated from flight data at different values of the angle of 
attack between 10 and 50 . The main part of the data was obtained from 
small amplitude longitudinal and lateral maneuvers. A small number of large 
amplitude maneuvers was also used in the estimation. The measured data were 
first checked for their compatibility. It was found that the accuracy of 
air data was degraded by unexplained bias errors. Then, the data were 
analyzed by a stepwise regression method for obtaining a structure of 
aerodynamic model equations and least squares parameter estimates. Because 
of high data collinearity in several maneuvers, some of the longitudinal and 
all lateral maneuvers were reanalyzed by using two biased estimation 
techniques, the principal components regression and mixed estimation. The 
estimated parameters in the form of stability and control derivatives, and 
aerodynamic coefficients were plotted against the angle of attack and 
compared with the wind tunnel measurements. The influential parameters are, 
in general, estimated with acceptable accuracy and most of them are in 
agreement with wind tunnel results. The simulated responses of the aircraft 
showed good prediction capabilities of the resulting model. 
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vector of input variables 
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j-th eigenvalue of X X matrix 

vector of residuals 
3 

air density, kg/m 
correlation coefficient 


c.g. 
m. a. c. 
Subscript: 
E 

ME 

PC 


center of gravity 
mean aerodynamic chord 

measured value 
mixed estimator 

principal components estimator 
initial value 


Superscript: 


unbiased estimated value 
biased estimated value 
Matrix exponent: 


T 

-1 


transpose matrix 
inverse matrix 


iv 


Derivatives of aerodynamic coefficients CMa 


Y, Z, l t m, n) 
referenced 

system of wind axes with the origin at the airplane center of gravity: 


referenced to a system of body axes and derivatives of referenced to a 


C 

a 

P 


a 


® j 


dc 

a 

p — 

dC 

a 

ac 

r - a 

oPb 

% 

L/ " -p 

a 

q 

d 2V 

a _ fl rb 

r a 2V 

dC 

a 

p — 

dC 

a 


da~ 

a 3 

W 


dc 

a 

! *V 

j = A, a, 

dh, f. 

h, £f , r 


v 




INTRODUCTION 


NASA began the High Alpha Technology Program (HATP) in 1988 with the 
main goal to accelerate the development of technologies which would expand 
high-angle-of -attack capabilities of future fighter aircraft. The flight 
research portion of the program has been using the High Angle-of-Attack 
Research Vehicle (HARV) which is a modified F/A-18 aircraft. One of the 
objectives of the flight program is to obtain low-speed, 
high-angle-of-attack aerodynamic parameters. These parameters may be used 
for validation of theoretical and wind tunnel predictions and for the 
development of a mathematical model of aircraft aerodynamics. 

As reported in references 1 and 2, the first attempts to obtain 
accurate parameter estimates of the HARV from flight data were not very 
successful. The main reason was the insufficient excitation of response 
variables in maneuvers intended for parameter estimation. Low excitation of 
these variables was caused by poor selection of input forms and, in some 
cases, by difficulties in maneuvering the aircraft in the requested way. 
The accuracy of estimated parameters was further degraded by close 
relationships between deflections of various control surfaces introduced by 
the HARV control system, uncertainty in the model structure at high 
angles of attack, and inaccuracy of measured air data. 

To avoid some of these problems, a new set of inputs was selected and 
verified in the flight simulator. Increased attention was given to the 
accuracy of measured incidence angles and to the extensive use of system 
identification methodology. 

The purpose of this report is to present a quick release of parameter 
estimates from various sets of data, compare those estimates with the 
existing aerodynamic model in the flight simulator, and demonstrate the 
prediction capabilities of the model determined from flight data. The 
report starts with the description of the aircraft, and flight and wind 
tunnel data available. Then, procedures for data analysis are briefly 
outlined. The results that follow include checks on the compatibility of 
measured responses, variation of estimated parameters with the angle of 
attack, and comparison of these estimates with wind tunnel measurements. 
The existing inconsistencies in flight results and discrepancies between 
them and wind tunnel data are discussed. Finally, measured and predicted 
aircraft motion in small amplitude maneuvers are compared. 
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AIRCRAFT 


The test vehicle is a twin-engine, single-seat fighter aircraft. It 
has a moderately swept wing with highly swept leading-edge extensions. The 
all-moving horizontal tail surfaces are mounted behind and below the wing; 
twin vertical tails are canted and toed out. The aircraft is controlled by 
four digital computers working in parallel. The computers are used in 
conjunction with redundant electrohydraulic servoactuators and analog 
sensors to provide primary control capabilities. There is also a backup 
mechanical control of the stabilator surfaces and open-loop analog control 
of the aileron and rudder. Longitudinal control uses symmetric deflections 
of the stabilator, leading edge and trailing edge flaps. Lateral control is 
provided by the ailerons, differential deflections of the stabilator, 
leading and trailing edge flaps, and synchronous rudder deflection. A 
drawing of the aircraft is presented in figure 1. The basic geometric, 
mass, and inertia characteristics are summarized in table I. A more 
detailed description of the aircraft and its control system is contained in 
references 3 and 4. 

The test aircraft was modified by adding right- and left-wing-tip booms 
with Pitot-static heads and a- and (3-vanes as shown in figure 1. 
Furthermore, flush pressure orifices were mounted on the forward radome area 
of the aircraft and on the remaining part of the forebody and leading-edge 
extension. The air data could be obtained from pressure measurements on the 
radome. The aircraft has a pulse-code modulation instrumentation system 
with telemetry as the only source of data. The measured data are recorded 
at the telemetry ground station. The instrumentation system includes 
transducers for the measurement of closed- and open-loop input variables, 
response variables, control system and engine operation, and fuel 
consumption from which instantaneous mass and inertia characteristics were 
calculated. 
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FLIGHT AND WIND TUNNEL DATA 


The flight data of the tested aircraft were obtained from NASA Dryden 
Flight Research Facility in the form of time histories sampled at 50 
samples/sec. The measured data were corrected for the c.g. offset of the 
linear accelerometers and wind vanes, and for the upwash and sidewash 
effects of the a- and 0-vanes. The air data for the analysis were taken as 
the average values from the right- and left-boom sensors. Various maneuvers 
were initiated from mostly steady flights at altitudes between 5,000m and 
12,000m (17,000 ft and 39,000 ft). In all maneuvers the Mach number did not 
exceed the value of 0.4. The scheduling of leading and trailing edge flaps 
with the angle of attack is shown in figure 2. 

Three different sets of maneuvers were available for the analysis. The 
first set consisted of 32 longitudinal and 32 lateral small amplitude 
maneuvers at angles of attack between 10 and 50 . The pilot input for the 
longitudinal maneuvers was a pitch command usually applied as three doublets 
of various duration. For the lateral responses, separate yaw and roll 
commands in the form of simple doublets were used. Time histories of 
open-loop input and response variables from the maneuvers are presented in 
figures 3 and 4. In figure 3 both the longitudinal and lateral maneuvers at 
a * 13° are shown. The inputs for the longitudinal maneuver included 
deflections of the horizontal tail, leading edge flaps, and trailing edge 
flaps. In the lateral maneuvers the input variables were the aileron, 
rudder, and differential tail deflections. Because of low airspeed (below 
120m/sec) the control system did not move the differential leading and 
trailing edge flaps. Both maneuvers represent good excitation of all 
response variables. As expected, the control system introduced strong 
coupling between deflections of aileron and differential tail, and coupling 
between symmetric leading and trailing edge flaps deflection, and angle of 

O 

attack. Figure 4 shows an example of similar maneuvers at a « 44 . From 
the time histories, problems of insufficient excitation of linear 
accelerations and maintaining uncoupled responses were visible. Because of 
the high values of a in these maneuvers the leading and trailing edge flaps 
remained in fixed position. On the other hand, strong coupling introduced 
by the control system existed between the aileron, differential tail, and 
rudder . 

The second set of data included 14 large amplitude maneuvers; 4 
longitudinal and 10 lateral. In each of these maneuvers, the motion^was 
excited within the extended range of angle of attack, usually from 10 to 
50°, using commanded doublets of various amplitudes and durations combined 
with a gradual Increase of the horizontal tail deflection. One of the large 
amplitude lateral maneuvers is shown in figure 5. Finally, the last data 
set contained three quasi-steady deceleration-acceleration maneuvers. 

The wind tunnel data of the F/A-18 aircraft are summarized in 
reference 3. The stability and control derivatives denoted as "wind tunnel" 
in this report were computed from aerodynamic functions used in the NASA 
Langley Research Center flight simulator of HARV under the following 
conditions: the Mach number of 0.4, altitude of 20,000 ft, c.g. location at 
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25 percent of the m.a.c., scheduled flaps position and horizontal tail 
deflection required to trim the aircraft at given angle of attack. The 
angles of attack varied between 2 to 54 with the increments of 2 °. The 
aerodynamic coefficients C L> C D and C m were obtained as functions of a for 

the above mentioned conditions but referred to 6 = -6°. 


FLIGHT DATA ANALYSIS 

The first step in the analysis included a check on measured data 
compatibility and estimation of unknown bias errors in the measurements. 
Then, a structure of aerodynamic model equations was determined and unknown 
parameters, mostly in terms of stability and control derivatives, estimated. 
The accuracy of least squares estimates can be, however, degraded by 
near-linear dependency (collinearity ) among measured time histories. The 
existence of data collinearity and its possible effect on the estimates 
were, therefore, estimated. As the result of that, some maneuvers had to be 
analyzed again using different estimation techniques which can reduce 
damaging effect of collinearity. 

For the compatibility check the maximum likelihood method of 
reference 6 was applied. The state equations were represented by kinematic 
equations 

x = f ( x , T) , 0 ) (!) 

where 

x = [u,v,w,^ ( e,0] T 
~ t a x . a y .a z >P.q.r] T 

and 0 is a vector of unknown bias and scale factor errors in measured input 
and response variables. The vector of response variables was formulated as 

z = [V,£,a,#,0 ,(/»] T 

and each measured response variable was expressed as 

z =(l+A)z + b+ e (2) 

t z z z 

where A z is the unknown scale factor error, b is the constant bias error 
a nd e z is the measurement noise. For the measured inputs 17 it was assumed 
that the scale factor errors and the measurement noise are equal to zero. 
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The unknown parameters and their Cramer-Rao bounds were obtained by 
minimizing 


N 

J = - i £ v T (i) R _1 v(i) - " n | R | (3) 

i=l 


where 


i>(i) = z^U) " z(i,0) 

R is the covariance matrix of measurement noise and N is the number of data 
points. 

A stepwise regression of reference 7 was used for model structure 
determination and parameter estimation. In linear regression the 
aerodynamic model equations are formulated as 

y = 0 + 0 .X. + ... 0 X ( 4 ) 

J 0 1 1 n n 


where y is the aerodynamic coefficient (dependent variable), x^ to x^ are 

the measured response and input variables or their combinations (regressors) 
and Qy j = 0, 1, . 


. . , n are the unknown parameters. The aerodynamic 

coefficients were calculated from the following expressions 
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In these equations the thrust was computed from the engine subroutine in the 
F/A-18 simulator for given M, h, and power lever angle. The angular 
accelerations were obtained by fitting cubic splines to measured angular 
velocities and subsequent differentiation of fitted curves. 

Candidate regressors for the small amplitude maneuvers were postulated 
as 

a, qc/2V , 5 h> 5 f , 6 ff , a 2 , a 3 , aqc/2V 

aS h , 0 2 , V/V Q , | rb/2V | , |pb/2V| 

and for the small amplitude lateral maneuvers as 

0. pb/2V , rb/2V, 5 , 6 , 0 2 , 0 3 

a r 

afi, apb/2V, arb/2V, aS , a<5 

a r 

<x 2 0, a 2 pb/2V , a 2 rb/2V, a 2 5 , a 2 8 , a 

a r 

Because of linear relationship between 5 and 5 „ , the differential tail 

a dh 

deflection was combined with aileron deflection to introduce the following 
control effectiveness term 




a = V, l, or n 


where 5^/5^ was estimated from measured data as 0.420 ± 0.0066. 


(6) 


The large amplitude maneuvers were analyzed after partitioning an 
ensemble of data from repeated measurements into subsets of selected 
a-intervals, as described in reference 8. The distribution of data points 
in these subsets for longitudinal and lateral maneuvers is shown in 
figure 6. The regressors for the longitudinal and lateral subsets were 
postulated as 


qc/2V, 3 h , 3 f , 3 <f . 3 2 . 3 2 , 6= 

and 

0, pb/2V , rb/2V , 5 6 0 3 , 06. , (pb/2V) 2 , (pb/2V) 3 

a r n 


respectively. 
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Finally, for a single large amplitude longitudinal maneuver a model was 
formulated as 


C = C (a), ,o + c (a) §§ + C (a) A5, + C 
a a 5 = -6 a 2V a s . h a 0 2 

h q 6h 0 


(7) 


where 


Ad, = 5 - 


-6 


h h 57.3 h 


= + 0.1047 


The first three terms on the right hand side of eq. (7) were approximated by 
the first-order polynomial splines (see reference 9). The control terms 

and A5^ were not included because of their small effect on the estimates of 

remaining parameters. 


The unknown parameters in eq. 
function 



(4) were obtained by minimizing the cost 


x j ai 


4 


(8) 


where l is the number of statistically significant terms in eq. (4). The 
least squares estimates of unknown parameters were obtained as 

0 = (X T X) X T Y (9) 

where X is the matrix of regressors and ones, and Y is the vector of 
measured dependent variables. The covariance matrix of parameters was 
estimated as 

Cov (0) = s 2 (X T X) (10) 


2 

where s is the variance of the measurement noise. The square root of the 
variance can be interpreted as a fit error for aerodynamic coefficients. 
Its estimate is based on the residuals and has the form 


j 


N 

E ty F (i) - y<n]‘ 

i=i b 

FI 


( 10 ) 


In some cases t-statistics (the t-distribution) should be included if the 
interest is focused on the significance of parameters. The estimates of 
these statistics are given as 


* 6 J 

” sTO 


(ID 
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The possibility of data col linearity in measured data was investigated 
by procedures described in reference 10. They include 

T 

a) Examination of the correlation matrix X X where the regressors are 
standardized (centered and scaled to unit length). 

T 

b) Eigensystem analysis of the X X matrix. The eigenvalues close to 
zero indicate near-linear dependency in the data. As a measure of 

T 

the spread of the eigenvalues of X X, the condition number, 
defined as the ratio of the maximum to minimum eigenvalue, is 
used. Condition numbers between 100 to 1000 imply moderate to 
strong collinearity. 

c) Parameter variance decomposition into a sum of components, each 

T 

corresponding to one, and only one, of the eigenvalues of XX. An 
unusually high proportion in the variance of two or more 
parameters for the same small eigenvalue can provide evidence that 
the near dependency is causing problems. 

The application of the ordinary least-squares technique to a set of 
collinear data very often results in nonphysical values for parameters and 
large values of their covariance. In order to obtain more stable and 

accurate estimates, two biased estimation techniques of reference 10, the 
principal components regression and mixed estimation, were applied to 
maneuvers where data collinearity was detected. These techniques provide 
estimates which are biased but have smaller variance than that of the least 
squares estimates. 

The principal components regression technique uses orthogonal 
regressors rather than the original ones. The orthogonal regressors are 

arranged in order of decreasing eigenvalues of X^X. Then, the last n of 

r 

these eigenvalues are removed from the analysis and the least squares 
principle is applied to the remaining components. Then, the parameter 
estimates associated with the orthogonal regressors are transformed back to 
the set of original parameters. The principal components estimator of 0 
thus takes the form 


n+l-n 


0 PC = E r i- t* x T vt 


j=i 


a j J 


■j 


and the covariance matrix has the form 

n+l-n 


cov«e pc ) - s 2 £ r i- t jt / 


(12) 


(13) 
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where Aj are the eigenvalues of the X X and t^ is the eigenvector 
corresponding to A. 

The assumption of an integral rank for X can sometimes be too 
restrictive, especially if the number of regressors is small. A possible 
improvement to the principal components estimator, known as the fractional 
rank estimator, has been proposed in reference 11. 

The mixed estimation is a procedure which uses prior information 

d * P0 + C (14) 

to augment the measured data. The mixed estimator is obtained as 


- [Ij x T x * pV'p]" 1 [1j x T v . pV'd] 


with covariance matrix 


cov(9 me ) 


IT T - 
XX + F W 

-s 


w-vr 


In eq.(14) to eq. (16) ( is a vector of random variables with E(£) * 0 and 
T 

E(CC J = W, and P is a matrix of known constants. 
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RESULTS AND DISCUSSION 


The results of the flight data analysis are summarized in the following 
six sections. The first one describes the estimates of bias errors in 
measured data and the compatibility between measured and predicted time 
histories of response variables. In the second section, structures of 
adequate models for longitudinal and lateral maneuvers are presented. This 
section is followed by investigation of data collinearity and by measures 
taken to avoid its damaging effect on parameter estimates. The next two 
sections include estimated longitudinal and lateral parameters and their 
comparison with wind tunnel data. The last section presents three examples 
demonstrating prediction capabilities of the resulting model. 

Data Compatibility Check: 

A selected number of longitudinal and all lateral maneuvers were 
analyzed to obtain estimates of bias errors in the measured response 
variables, and residuals represented by differences between measured and 
predicted response variables. The parameter estimates indicated a need for 
additional corrections to the measured data, and the residuals provided 
information about bounds on the remaining errors in the measured data. The 
results from the longitudinal and lateral maneuvers in figure 3 are 
presented as an example. The parameter estimates and their standard errors 
are given in tables II and III; the time histories of output variables and 
corresponding residuals in figures 7 and 8. The estimates indicate small 
errors in angular rates, Euler angles and angle of attack, larger errors in 
longitudinal and lateral accelerations, and an excessive scale factor error 
in the sideslip angle. Similar conclusions could be drawn from the 
remaining maneuvers. Because the analysis of small amplitude maneuvers uses 
the increments of measured variables, the main concern was the large scale 

factor error in sideslip angle. The estimates of A from various maneuvers 
( P 

varied randomly with no apparent dependence on the angle of attack or 
magnitude of sideslip angle. The average value of A„ obtained from 32 small 

R 

amplitude maneuvers was equal to 0.130 with an ensemble standard error of 
0.097. Each sideslip angle time history was corrected by using the estimates 
corresponding to that maneuver. The estimates of the scale factor 

errors for a also varied randomly around zero not exceeding the value of 
0.02. For that reason no corrections to the angle of attack were applied. 

The residuals in figures 7 and 8 still include the effect of 
uncorrected bias errors. The magnitude of these errors is mainly visible in 
the air data variables. By examining all maneuvers it was found that the 
bias errors in measured air data have bounds equal to ±3 m/sec for the 
airspeed and ± 1 for both a and 0. In some maneuvers for a > 40° the 
bounds for V and 0 were increased to ± 5 m/sec and ± 2.5°, respectively. 
The bias errors in the air data are usually related to one short segment of 
the whole maneuver. Despite this, the quoted inaccuracy in the air data 
measurement is unacceptable for a serious research test. The main factor 
contributing to these errors is the location of the vanes and Pitot-static 
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heads. Because this arrangement cannot be changed, it is necessary in 
future data compatibility checks to consider changes in the postulated 
model. Possible changes may be to include measured position of the aircraft 
and/or to include more unknown parameters in the measurement equation. The 
other possibility is to use the flush air data system after its performance 
is checked. 


Adequate Models for Aerodynamic Coefficients: 

A stepwise regression method was applied to each set of transient data. 
As a result, adequate models for the aerodynamic model equations were 
determined and the least squares estimates of parameters in the model 
obtained. The selection of these models was based mainly on changes in the 

2 

multiple correlation coefficient, R , with an increasing number of terms 
included in the model. In addition to this criterion, the fit error and 
statistical significance of the estimates were examined and the residuals 
checked for unexplained differences between measured and predicted 
aerodynamic coefficients. More about the selection of an adequate model can 
be found in reference 7. An example of measured and predicted aerodynamic 
coefficient is given in figure 9. These coefficients correspond to the 
lateral maneuver shown in figure 3. In this example the residuals reflect 
the effect of rather excessive measurement noise and model inadequacy. 


A general form of an adequate model for longitudinal small amplitude 
maneuvers was determined as 


C 

a 



2 3 

+ C Aa + C 0 Aa + C - Aa 
a a 2 a 3 

a a a 


+ 



a 

q 




+ c 


AS, 


Sf 


(17) 


+ C AS fljr , a = L, Z, or m 

a 5if tf 

The reference conditions in this model correspond to trimmed a, scheduled 
flaps position and 5^ = “6 . Therefore, 


Aa = a - <Xq 



A5 £f & If 8 lf Q 


AS, = S, + 0.1047 
n h 
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In many maneuvers, especially those with small excitation, a linear model 
was found to be the best. The leading and trailing edge flap derivatives 
appeared significant only in a limited number of maneuvers for a between 10° 
to 26 . 


Adequate models for lateral small amplitude maneuvers have the form 
C 


C +C p + C„3 3 + C E* + C — 
a„ a a 3 a 2V a 2V 

P P P r 


‘0 


+ C 8 + C 8 

a 5 a a 5 r 
a r 


a = Y, t, or n 


(18) 


The cubic term in £ was significant only in models for the rolling-moment 
coefficient with a between 35 to 40°. 


Models for partitioned data usually Included linear terms only. Many 
stability and control derivatives could not be estimated because of low 
significance of associated regressors. For a single large amplitude 
longitudinal maneuver, first-order splines were adequate for the 
approximation of coefficients and derivatives. As in the case of 
partitioned data, low excitation of the short-period mode prevented more 
detailed determination of model structure. 

Effect of Data Collinearity on Parameter Estimates: 

After determining model structures and estimating parameters, several 
small amplitude maneuvers were checked for data collinearity and its 
possible effect on parameter estimates. Two examples of these checks are 
given. The first one uses data from the longitudinal maneuver in figure 3. 
An adequate model for the pitching-moment coefficient includes eight 

regressors. The X X matrix in correlation form and eigenvalues of the X^X 
matrix are presented in table IV. Two high pairwise correlations (a 8 ) 

3 ’ £f 

and (a, a ) exist. The condition number of 647 indicates a moderately 

strong effect of data collinearity on parameter estimates. For that reason 
the covariance decomposition proportions were not included in table IV. The 
least squares parameter estimates are summarized in table V. They are 
compared with the results of mixed estimation using C = 0.083 ± 0.013 as 

m S£f 

a loose a priori value. Table V also includes the standard errors and 
t-statistics of parameters and fit errors, sfC^), for both techniques. From 

the comparison of results it follows that all parameters in the model, with 
the exception of the least squares estimates of C , are statistically 

m 6£f 

significant (t* > 1.96). The differences between the two sets of parameters 
are within 2<r-bounds, but the introduction of a priori information sharpens 
the regression results. Based on this observation the maneuvers with 
leading- and trai ling-edge flaps active were analyzed again using the mixed 
estimation technique. The new parameter estimates represent the final 
results discussed in the next section. 
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In the second example, the data collinearity diagnostic for a lateral 

maneuver at a « 41° is presented in table VI. In this maneuver strong 

correlation exists between the aileron and rudder deflection. The spread of 

eigenvalues is very small, however, resulting in a condition number equal to 

47. The effect of strong pairwise correlation p(5 , 6 ) on parameters in 

a r 

the rolling-moment equation is demonstrated in table VII by using three 
different estimation techniques. The results from consecutive parameter 
entries of the stepwise regression algorithm show the change in the value of 
C, from -0.041 to -.102 when the regressor 5 enters the model. This 

sudden change is a direct result of strong correlation between 8^ and 6^. 

When the mixed estimation with an a priori value of C. = 0.000 ± 0.010 was 

Sr 

introduced, the value of C. = 0.055 was obtained. This value agrees well 

oA 

with the wind tunnel prediction of C. . The estimates of C. were further 

SA oA 

verified by using the principal components regression. The reduction of 
T 

rank of matrix X X from 6 to 5 was too coarse, as shown by the resulting fit 
error. Therefore, several partial rank reductions were tried. For a rank 
of 5.5, both the fit error and parameters were almost identical to those 
from the mixed estimation. The t-statistics indicate that, in all cases, 
the parameters and have low significance in the model. In this 

P r 

example an adequate model would have only three regressors, /3, 8 , and 8 , 

a. r 

which is also indicated by the values of the squared multiple correlation 
coefficient. 

Because of the damaging effect of the aileron-rudder correlation on 

parameter estimates, the collinearity diagnostic was performed on all 

lateral maneuvers. The resulting correlation coefficient, p(8 , 8 ), is 

a r 

plotted against the angle of attack in figure 10. It can be expected that 
the least squares estimates of the control parameters from maneuvers at a > 
20° can be influenced by this correlation. All the lateral maneuvers were, 
therefore, analyzed again and the mixed estimates were selected as the final 
set of parameters. The a priori values of C y , CL and C were based on 

SA or n 8A 

wind tunnel data. The uncertainty in these values was expressed by standard 
errors varying between 0.02 to 0.04. The level of uncertainty was kept as 
high as possible depending on the correlation between estimated control 
parameters. Values of this correlation coefficient less than 0.85 were 
considered acceptable. 

Longitudinal Parameters: 

The longitudinal parameters include three coefficients, C^, and C^, 

and their derivatives plotted against the angle of attack. The three 
coefficients corresponding to steady conditions with scheduled flaps and 
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fixed horizontal tail deflection, 5^ = -6 , are presented in figures 11 to 

13. In figure 11 the wind tunnel measurement of is compared with flight 

results from a slow deceleration/acceleration maneuver. The values from 
both parts of the maneuver agree quite well, but for a between 15° to 45° 
they are slightly lower than those from the wind tunnel measurement. The 
flight data of and in figure 12 contain the results from small 

amplitude maneuvers and from a single large amplitude maneuver. The 

agreement between both sets is very good. The flight data of C L are closer 

to the wind tunnel curve than the values of C z which have lower values 

almost over the whole range of a. Figure 13 presents the C coefficients 

m 

estimated from small amplitude maneuvers, partitioned data and single large 

amplitude maneuvers. All three sets of flight results agree well. Some 

differences exist, however, between flight and wind tunnel measurement in 

the region of a between 20° and 40°. 

The derivatives of longitudinal coefficients with respect to a are 
plotted in figure 14. They agree, in general, with the wind tunnel 

predictions. The o only pronounced disagreement can be seen in the parameter 
around a = 40 where the flight data do not indicate a sudden increase 
a 

of static stability as the wind tunnel test does. The q— derivatives are 
given in figure 15. Both the C L and C z parameters are in sharp 

q q 

disagreement with the wind tunnel data for a < 40°. These differences have 
only a small effect on the prediction of aircraft motion because of the low 
significance of the corresponding terms in the model. Parameter values of 
C m from sma11 amplitude maneuvers and partitioned data indicate a sudden 

q 

increase in damping at a ~ 22° which was not predicted by wind tunnel 
measurement. Different values of C from the large amplitude maneuver for 

o q 

a between 12 to 22 might be the result of limited information in the data 

which prevented more refined determination of C (a) variation. For 

m 

q 

apparently the same reason no estimates of C for a > 38 were obtained 

m 

q 

from partitioned data. 


Figure 16 presents control derivatives related to horizontal tail 
deflection. All three parameters indicate lower tail effectiveness than 
that predicted from wind tunnel measurement. Large inconsistency exists in 
the C m estimates from three different sets of data. The differences are 
oh 

sometimes greater than 0.1 which is greater than 15 percent of the estimated 
values. The large scatter in the estimates is surprising because the term 
c m A6 is the most significant in the model which means that C should be 
Sh n m 5h 
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very well identifiable. At the same time the standard errors of the C 

m Sh 

estimates were low, varying between 0.005 to 0.014. No reasons for the low 
accuracy of C parameters were found. 
m 3h 

The parameters expressing the effect of flaps are given in figures 17 
and 18. The trai ling-edge derivatives confirmed the trend in these 
parameters with the angle-of-attack change. The closeness of the 
leading-edge derivatives to the wind tunnel data is caused by their 
selection as a priori values in the mixed estimation and, at the same time, 
their low sensitivity. 

Lateral Parameters: 

The lateral stability and control parameters obtained from flight and 
wind tunnel measurements are included in figures 19 to 23. All the lateral 
parameters in a linear model were estimated from small amplitude maneuvers 
with different levels of consistency. The scatter in the estimates 
indirectly indicated the significance of corresponding terms in the model 
equations. Only the most important parameters were estimated from 

partitioned data, because the remaining parameters did not enter the model 
during the stepwise regression analysis. 

The flight and wind tunnel sideslip derivatives are compared in 
figure 19. The differences in C for a less than 35 and greater than 40 , 

p 

and the differences in C. around a = 22 could not be explained. The 

p 

decrease in directional stability in flight for a > 40 is caused by the 
horizontal tail setting. The wind tunnel data correspond to 6^ = 0 whereas 

o 

during the flight maneuvers at a > 40 the tail deflection varied around 
5 h = - n °- 


Figure 20 contains the derivatives with respect to roll rate. The 

P 

parameter varies substantially around the wind tunnel data with a tendency 
for large negative values for a > 40 . The parameter exhibits gradual 

P 

decrease in the roll-damping with the increase of a. The difference between 

flight and wind tunnel results in the region of a between 25 to 42 could 

be the result of different types of maneuvers. Flight data were obtained 

from transient maneuvers while the wind tunnel testing was based on forced 

oscillations with small amplitude in roll angle equal to 5 . The variation 

of cross-derivative C for a less than 38 is similar to that from the wind 

n 

P 

tunnel test. For a > 40°, however, the flight data indicate an increase in 
negative values of this parameter. 
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The yaw-rate parameters are included in figure 21. Low excitation of 
yawing velocity in transient maneuvers and low contribution of yaw-rate 
terms to the lateral motion is reflected in large scatter of all three 
parameters. The resulting parameters are close to the wind tunnel 
prediction for a < 40 . For a > 40° there is a change in Cy to higher 

r 

values and in from positive to negative value. All the control 

r 

parameters in figures 22 and 23, except C , were estimated with high 

n 6A 

consistency and closeness to the wind tunnel data. The values of the 
adverse-aileron-effect parameter, C , are consistently lower than those 

n 6A 

from wind tunnel measurements. 

Model Prediction Capabilities: 

The prediction capabilities of the model determined from flight data 
were checked in several longitudinal and lateral maneuvers not used in the 
previous analysis. For each type of motion, the simulation was based on 
models with three degree-of-f reedom. The remaining variables in the 
equations of motion were substituted by measured values. 

One of the simulated longitudinal maneuvers at a « 32° is presented in 
figure 24 where the time histories of the input variable and three output 
variables are plotted. The frequency of the motion is predicted very well. 
There are some discrepancies in the amplitudes of all three output 
variables. The remaining two figures include the predicted and measured 
lateral maneuvers. As in the previous case, the time histories of input and 
output variables are given. The first maneuver presented in figure 25 
compares the measured and predicted output variables excited at a « 12°. 
The agreement between the two sets of data is very good. The second 
maneuver in figure 26 was performed at a ~ 48 . In this high-a maneuver the 
prediction differs from the measurement, but the main features of the 
response are predicted quite well. The deterioration of model prediction 
capabilities with increasing angle of attack was also observed in the 
remaining simulated maneuvers. 
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CONCLUDING REMARKS 


Aerodynamic parameters of the High-Angle-of-Attack Research Vehicle 
(HARV) were estimated from different types of maneuvers at angles of attack 
varying between 10 and 50 . Data analysis included the data compatibility 
check, model structure determination, and parameter estimation. For these 
three steps a maximum likelihood estimation, stepwise regression, and three 
additional estimation techniques: the ordinary least squares, principal 

components regression, and mixed estimation, were applied. The resulting 
parameter estimates were obtained in the form of stability and control 
derivatives and aerodynamic coefficients. They were then presented as 
variations with the angle of attack and compared with wind tunnel 
measurements. In addition to the data analysis, the model prediction 
capability was checked. From all the results obtained the following 
conclusions can be drawn: 

1. The data compatibility check showed good agreement between measured and 
predicted attitude angles. Some unexplained errors existed, however, 

in the air data. The error bounds for a and 0 were ± 1° with an 

o o 

increase up to ± 2.5 for 0 in some maneuvers with a above 40 . The 
effect of these errors on parameter estimation was not investigated. 

2. Most of the adequate models for aircraft aerodynamic coefficients 
contained only linear terms with stability and control derivatives. 
For the longitudinal maneuvers with large excursions and lateral 

o o 

maneuvers for angles of attack between 25 and 42 , the adequate models 
also included nonlinear terms in a and 0. 

3. Increased data collinearity was present in longitudinal maneuvers with 
extensive excitation and leading edge and trailing edge flaps 
operating. In lateral maneuvers the near linear dependency existed 
among the aileron, differential tail, and rudder deflections. For the 
data with high collinearity, the mixed estimates of parameters were 
found more accurate than the estimates obtained by the ordinary least 
squares technique. 

4. The lift, vertical-force, and pitching-moment coefficients estimated 
from different types of maneuvers as functions of angle of attack are 
in general agreement with wind tunnel tests. 

5. From the measured data it was possible to estimate all stability and 
control derivatives. However, the accuracy of these estimates and 
their closeness to wind tunnel data vary substantially. Lateral 
control parameters are estimated with high consistency and are very 
close to the wind tunnel data. Larger scatter than expected existed in 
the estimates of the horizontal tail effectiveness. These parameters 
also have lower values than those predicted by wind tunnel measurement. 
All static stability parameters are estimated with acceptable accuracy 
expressed by their scatter. Unexplained departures from wind tunnel 
data exist in the longitudinal static stability parameter for angle of 
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o 0 

attack around 40 , dihedral-effect parameter around a = 22 , and in the 
lateral force parameter for almost the entire range of angle of attack. 
Almost all the dynamic parameters exhibit large scatter due to their 
low sensitivity. The most important differences between wind tunnel 
and flight data are in the damping-in-roll parameter and 
rolling-moment-due-to-yawing parameter. The flight results do not 

indicate an increase in roll damping between angles of attack of 25° 

O 

and 42 as shown by wind tunnel data. Values of rolling moment with 
yawing velocity, , change from positive to negative for angles of 

r 

attack greater than 40°. 

The future experiment and data analysis should include the following: 

1. Maneuver for the assessment of data compatibility and accuracy using 
air data from the wing-boom sensors and flush air data system. 

2. Different types of maneuvers for obtaining parameters with increased 
accuracy. 

3. Evaluation of measured data accuracy and different estimation 
techniques on the accuracy of parameter estimates. 

4. More theoretical development in techniques for model structure 
determination and parameter estimation from flight data with high 
collinearity. 
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Table I. Geometric, Mass and Inertia Characteristics 
of the Aircraft 


Total Length, m 

17.07 

Wing: 2 

Area, m 

37.16 

Span, m 

1 1.41 

Mean geometric chord, m 

3.51 

Aspect ratio 

3.5 

Quarter-chord sweep angle, deg 

20.0 

Horizontal tail: 2 

Area (wetted), m 

16.35 

Span, m 

6.58 

Mean geometric chord, m 

1.91 

Aspect ratio 

2.40 

Quarter-chord sweep angle, deg 

42.8 

Moment arm (e g. at 0.25 m.a.c.), m 

5.12 

Vertical tail: 2 

Area (wetted), m 

9.66 

Mean geometric chord, m 

2.13 

Aspect ratio 

1.20 

Quarter-chord sweep angle, deg 

35.0 

Cant, deg 

20.0 

Moment arm (e g. at 0.25 m.a.c.), m 

3.10 

Mass, kg 

14,400 

Inertia: 2 

x • kg™ 

28,880 

2 

ly t kg-m 

165,930 

2 

>Z * 

185,030 

2 

xz • k g-m 

-2,630 
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Table II. Estimates of Bias Errors in Measured 
Longitudinal Data 


Parameter 

0 

8(0) 

by, m/sec 

-.14 

.035 

X 

a 

-.010 

.0014 

b^, deg/sec 

-.0106 

. 00020 

b 0 . deg 

-.386 

.0039 


.0505 

00040 

b , g units 
ax 

. 1 52 

.0021 

b , g units 
az 

.061 

. 0023 


Table III. Estimates of Bias Errors in Measured 
Lateral Data 


Parameter 

0 

9(0) 


-0.54 

0.023 

V 

0.171 

0.0029 

bp. deg/sec 

-0. 1564 

0.00075 

deg/sec 

“0. 0729 

0.00071 

b ^. deg 

-0.087 

0.0094 


0.071 

0.0011 

X 

f 

0.063 

0.0013 

b , g units 
ay 

-0.111 

0.0027 
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Table IV. Collinearily Diagnostic for Longitudinal Maneuver 


X X matrix (scaled and centered regressors): 


a qc/2V 



_ *f 

•0 

2 

a 

3 

a 

1.000 . 223 

.016 

.229 

.775 

-.313 

.881 

1.000 

.003 

. 193 

-.206 

-.079 

.227 


1.000 

-.487 

.367 

. 575 

-.042 



1.000 

.007 

-.687 

.325 




1.000 

-.579 

.648 





1.000 

-.470 






1.000 


eigenvalues (scaled 

regressors) : 


3.190, 2.559, 

1.125, .588, 

.301 , 

.156, .0765, 

. 00493 
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Table V. Least Squares and Mixed Estimates 
Longitudinal Maneuver 



Least squares estimation 

Mixed estimation* 

Parameter 

A 

9 

3(0) 

t* 

0 

3(0) 

t* 

c 

m 

o 

.079 

.0016 

50.6 

.0852 

. 00096 

84.9 

C 

m 

a 

- .09 

.030 

3. 1 

- .06 

.017 

3.2 

C m 2 
ot 

- 1.90 

. 075 

25.4 

- 1.88 

.074 

25.5 

C m 3 
a 

-15. 

1.0 

14.9 

-14.6 

.99 

14.8 

C 

m 

q 

- 4.3 

.85 

5. 1 

- 5.4 

.49 

11.1 

c 

m «h 

- .816 

.0091 

89.5 

- .852 

.0070 

117.9 

c 

m 5f 

.099 

.0098 

10. 1 

. 107 

.0081 

13.2 

C 

m a«r 

- .04 

.024 

1.7 

- .07 

.013 

5.6 

s(C ) 
m 

- .00565 

- 

- 

. 00565 

- 

- 


•) a priori value C = - .083 ± .015 
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Table VI. Collinearity Diagnostic for Lateral Maneuver 


T 

X X matrix (scaled and centered regressors) 


0 

pb/2V 

rb/2V 

a 

a 

8 

r 

1.000 

- .083 

- .014 

- .009 

.058 


1.000 

.728 

.053 

- . 177 



1.000 

. Ill 

.045 




1.000 

.944 





1.000 


eigenvalues (scaled regressors): 

1.9888, 1.738, 1.172, .819, .239, .0426 


24 


Table VII. Least Squares and Biased Estimates, 
Lateral Maneuver 


Parameter 

St< 

epwise regression 

Mixed 

estimation* 

Principal 

components j 

CM 

H 

C 

n = 3 

n * 4 

n = 5 

n = 6 

r n a5 

r n * 55 

o 

o 

-.0022 

-.0005 

-.0000 

-.0000 

-.0000 

.0008 

.0010 

j .0008 

% 

-.103 

-.102 

-. 106 

i 

M 

o 

OH 

-.106 

103 

102 

[ 

-. 103 


[36.3] 

[42.9] 

[47.1] 

[47.3] 

[47.3] 

[45.2] 

[44.7] 

[44.8] 

c < 

- 

- 

- 

-.05 

-.03 

-.05 

-.09 

-.05 i 

p 













[3.0] 

[1.8] 

[2.1] 

[3.9] 

[1.9] | 

c t 

r 



- 

- 

-.05 

-.07 

-.05 

-.08 






[1.5] 

[1.5] 

[1.0] 

[1.6] 

c t 

6 A 

_ 

-.041 

-.102 

-.095 

-.095 

-.055 

-.019 

-.059 



[17-0] 

[44.2] 

[41.2] 

[41.2] 

[12.7] 

[14.0] 

[12.8] 

5r 

- 


.0145 

.0126 

.0128 

. 0033 

.0046 

.0047 




[27.9] 

[24.4] 

[28.8] 

[8.8] 

[17.2] 

j 

[4.6] 

S «V 

.00500 

. 00418 

.00396 

. 00394 

.00394 

. 00408 

.00429 

. 00413 

R 2 

1 

66.0 

76.2 

78.7 

78.9 

79.0 

- 

- 

- 


•) a priori value C. » .000 t .0010 
or 

Note: figures in brackets are t - statistics 



— 17.8 ft — »-| 



Figure 1. Three-view drawing of test aircraft. 
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Figure 2. Leading and trailing edge flap schedule 
with angle of attack. 
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Figure 3. Time histories of measured input and response 

variables in small amplitude longitudinal and lateral 
maneuver (a=13°). 
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Figure 3. Continued. 
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Figure 3. Concluded. 
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t, sec 


Figure 4. Time histories of measured input and response 

variables in small amplitude longitudinal and lateral 
maneuver (a=44°). 
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Figure 4. Concluded. 
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Figure 5. Continued. 
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Figure 5. Concluded. 
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o 10 20 30 40 50 


a, deg 


Lateral data 



0 10 20 30 40 50 


a, deg 

Figure 6. Number of data points in subsets using partitioning 
of data from large amplitude longitudinal and lateral 
maneuvers. 
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Figure 8. Concluded. 
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P(»a. 


Figure 



10. Variation of correlation between aileron and rudder 
deflection with the angle of attack. 
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Wind tunnel 

Flight: 

o deceleration 









0 20 40 60 

a, deg 


Figure 12. Comparison of lift and vertical-force coefficient at steady 
conditions estimated from flight and wind tunnel data. 
Scheduled flaps, 5^ = -6 . 
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Wind tunnel 

Flight: 

o small ampl. maneuvers 
partitioned data 



a, deg 


Figure 13. Comparison of pitching -moment coefficient at steady 
conditions estimated from flight and wind tunnel 
data. Scheduled flaps, 8h = -6°. 
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Wind tunnel 

Flight: 

o small ampl. maneuvers 
partitioned data 





Figure 16. Comparison of longitudinal control parameters estimated 
from flight and wind tunnel data. 
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Figure 18. Comparison of parameters expressing leading-edge 

flaps effect estimated from small amplitude maneuvers 
and wind tunnel measurement. 
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Wind tunnel 

Flight: 

o small ampl. maneuvers 
partitioned data 



<b 


O 


Wind tunnel 

o Flight 



0 20 40 60 

a, deg 

Figure 21. Comparison of yaw-rate parameters estimated from 
flight and wind tunnel data. 
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Figure 24. Comparison of measured longitudinal time 
histories with those computed. 
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Figure 25. Comparison of measured lateral time histories 
with those computed. a 0 * 12°. 
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Figure 25. Concluded. 




09 


o 8 * 7 “ °® P^nduioo asoqj qjiM 
sauojsiq auiij leiajei painseaui jo uosuBdui 03 9Z ainSuj 


oas ') 

21 8 o 

1 ■ i sr- 


s)jun 6 

<Xg 

r 

6 - 

9- 

e- 

n oas/6ep 
u ‘j 

E 

9 

6 

oe- 

02- 

oi- 

0 oas/6ap 

'd 

01 
02 
08 


* ' • ♦ * * » ^ 
\ . \ jP ♦ * *’ V* ; "'!'a ^vV. C( - 

-Xdy^ : » 





papipajd — J 21 

pajnseaw + 



Figure 26. Concluded. 
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